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ABSTRACT 

We present key results from deep spectra of the Diffuse Ionized Gas (DIG) halo of 
the edge-on galaxy NGC 5775. [NII]A6583 has been detected up to z « 13 kpc above 
the plane in one of two vertically oriented long slits - making this the spiral galaxy 
with the greatest spectroscopically detected halo extent in emission. Key diagnostic 
line ratios have been measured up to z ~ 8 kpc, allowing the source of ionization and 
physical state to be probed. Ionization by a dilute radiation field from massive stars in 
the disk can explain some of the line ratio behavior, but departures from this picture 
are clearly indicated, most strongly by the rise of [O III] /Ha with z. Velocities of the 
gas in both slits approach the systemic velocity of the galaxy at several kpc above the 
plane. We interpret this trend as a decrease in rotation velocity with with essentially 
no rotation at heights of several kpc. Such a trend was observed in the edge-on galaxy 
NGC 891, but here much more dramatically. This falloff is presumably due to the 
gravitational potential changing with z, but will also depend on the hydrodynamic 
nature of the disk-halo cycling of gas and projection effects. More detailed modeling 
of the ionization and kinematics of this and other edge-ons will be presented in future 
papers. 



Subject headings: galaxies: individual (NGC 5775) - galaxies: ISM - galaxies: spiral 
galaxies: structure ~ stars: formation 



1. Introduction 

The Reynolds layer or Warm Ionized Medium of the Milky Way is the repository of about 
90% of the free electrons in the ISM. How this medium is ionized, how its large (~ 1 kpc; Haffner, 
Reynolds, & Tufte 1999) scale height is maintained, and how far it extends vertically are some 
of the main outstanding issues regarding this diffuse medium. New progress on these issues is 
being made with the Wisconsin Hq Mapper (WHAM; Reynolds et al. 1998). A complementary 



^Visiting Astronomer, National Optical Astronomy Observatories, Tucson, AZ 



- 2 - 



approach is to examine this medium in external galaxies, where it is generally referred to as Diffuse 
Ionized Gas (DIG). Narrow-band imaging, spectroscopy, and Fabry- Perot observations have, in 
recent years, allowed the brightness, spatial distribution, emission line ratios, and kinematics of 
DIG to be studied in many nearby galaxies (e.g. Golla, Dettmar, & Domgorgen 1994; Greenawalt, 
Walterbos, & Braun 1997; Hoopes, Walterbos, &; Rand 1999; Wang, Heckman, & Lehnert 1997; 
Rand 1997). 

Enhanced ratios of [SII]AA 6716,6731/Ha and [NII]AA 6548,6583/Ha in the Reynolds layer 
and external DIG layers relative to HII regions (e.g. Haffner et al. 1999; Golla et al. 1994; 
Ferguson, Wyse, k, Gallagher 1996; Rand 1997) can be explained if the dominant source of 
ionization is dilute radiation leaked out of star forming regions from massive stars. However, 
recent high-quality data indicate departures from this simple picture: first, the near constancy 
with z of [SII]/[NII] in both the Milky Way and external edge-on galaxies (Haffner et al. 1999; 
Golla et al. 1994; Rand 1998), whereas a rise is expected (e.g. Domgorgen & Mathis 1994; 
Sembach et al. 2000); second, [SIIJ/Ha and [Nil] /Ha reach values ^ 1 which are difficult for 
models to reproduce; and third, [OIII]/Ha in the edge-on spiral NGC 891 is found to rise with z 
(Rand 1998), whereas it is expected to fall as oxygen becomes predominantly singly ionized. 

Two explanations for this behavior have been put forth. While viable, both are limited by our 
knowledge of the energizing sources of halos in general. First, Reynolds, Haffner, & Tufte (1999) 
and Haffner et al. (1999) have pointed out for the Reynolds layer that the behavior of [SIIJ/Ha, 
[NII]/Hq; and [SII]/[NII] could be explained if gas temperature rather than dilution of the ionizing 
radiation field (and accompanying changes in the ionization state of the metals) were the key 
parameter that changed with z, since the first two ratios are very temperature sensitive while the 
third is not. This idea has yet to be applied to an external galaxy. The alternative explanation, 
which has generally been considered for external galaxies (e.g. Rand 1998; Galarza, Walterbos, & 
Braun 1999; Martin 1997), is a secondary source of ionization, such as shocks (e.g. Shull & McKee 
1979) or turbulent mixing layers (Slavin, Shull, &l Begelman 1993) which contribute a fraction 
of the Ha emission which is small but increases with z. The attractive feature of such ionizing 
sources is that they can yield a high [GUI] /Ha ratio without necessarily dominating the emission 
from the other observable lines. 

A second issue well addressed in external galaxies is the vertical variation of the DIG 
kinematics. It has been found in the DIG of NGC 891 that the observed velocities along one 
vertical slit through the disk become closer to the galaxy systemic velocity with increasing z, 
suggesting a fall in the rotation speed (Rand 1997). HI data suggest that the effect is widespread 
in the lower halo (Swaters, Sancisi, &: van der Hulst 1997). The rate of falloff will depend on the 
galactic potential and the hydrodynamic nature of the disk-halo gas cycle (Benjamin 2000). 

Finally, of great interest is the detectable extent of DIG layers. DIG in NGC 891 and NGC 
5775 has been detected to about z = 5 kpc (Rand 1997; Hoopes et al. 1999), and z = 6 kpc 
(Collins et al. 2000), respectively, while deep, wide-field imaging of NGC 4631 indicates emission 
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up to 16 kpc from the plane (Donahue, Aldering, & Stocke 1995). 

This paper presents spectra of the edge-on galaxy NGC 5775, which has been previously 
imaged in the Ho; line by Collins et al. (2000). It is an interacting galaxy (e.g. Irwin 1994) with 
a high far infrared luminosity and inferred far infrared surface brightness (Collins et al. 2000) 
indicating active star formation. Its DIG layer is very bright and extended, as mentioned above, 
and features some of the most prominent vertical filamentary structure found above a galactic 
disk. Here we focus on several dramatic new results from these spectra. A full analysis of these 
data and spectra of three other edge-ons in terms of the the two ionization/heating scenarios 
discussed above will be presented in a future paper (Paper II) by Collins & Rand (in preparation) . 
Another paper (Paper III) will examine the kinematics of the NGC 5775 halo. 

2. Observations 

The spectra were obtained at the KPNO 4-m telescope on 1999 June 10-13. The slit positions 
run perpendicular to the plane of the galaxy, as shown in Figure 1, where they are overlaid on the 
Ha image from Collins et al. (2000). The central positions of the slits are given in Table 1. The 
slit length is 5', the slit width is 2" for the galaxy observations, and the spatial scale is 0.69" per 
pixel. Slit 1 was chosen to include emission from the most prominent extraplanar DIG filament, 
while Slit 2 covers a region of weaker emission. The KPC-007 grating was used with the T2KB 
2048x2048 CCD, providing a dispersion of 1.42 A per pixel, a resolution of 3.5 A, and a useful 
coverage of 4700 — 68OOA. Many half-hour spectra were taken. Total integration times are given 
in Table 1. 

The basic reduction steps are as described in Rand (1998), and here we only describe 
processing of the final, stacked, calibrated spectra. Emission from some of the stronger lines 
extends nearly to the edge of the slit, complicating sky subtraction. In addition, there is a focus 
variation of unknown origin, causing the lines towards the ends of the slit to be significantly 
broader than those closer to the center, precluding accurate subtraction of the night sky lines. 
While the general sky background can still be corrected for, sky lines, when blended with galaxy 
emission lines, are deconvolved in the line fitting process. This blending was mainly a problem 
for the Ha line, which coincided with a blended pair of sky lines. Towards the ends of the slit, 
the combined intensity from these sky lines became very constant, indicating that there was no 
measurable Ha and that this intensity value could be used in determining the true Ha intensity 
elsewhere along the slit. A much fainter sky line near [O III] was corrected for in the same way. 
Strong sky lines also precluded measurements or meaningful upper limits on the [01] A6300 and 
HeIA5876 lines for Slit 1, and the [Nil] A5755 line for both slits. 

Line parameters reported here are determined from spectra averaged along the spatial 
direction, typically by 10 pixels, or 830 pc at an assumed distance of 24.8 Mpc (Irwin 1994). 
Line properties were determined with Gaussian fits, using a linear fit to the continuum on each 
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side of the line. For uncertainties on intensities and central wavelengths, the variance of these 
quantities along each slit for sky lines of a range of intensities were determined, and an estimate 
of the dependence of these variances on sky line intensity was formed. This relation was then used 
to determine uncertainties for the emission lines of interest. The detection limit for intensities is 
about 10~^^ erg cm~^ s~^ arcsec"^. 

3. Results 

One of the most important results is the tremendous height above the plane to which emission 
can be detected, providing unique information on the ionization state and kinematics of gas at 
large distances above the midplane. [N II] A6583 is detected to about z = 13 kpc and z = 7 kpc 
on both sides of Slits 1 and 2, respectively. Ha is detected up to about 2 = 8 — 9 kpc in Slit 1 
and z = 6 — 7 kpc in Slit 2. [O III] A5007 is detected to about z = 8 kpc on the NE side of Slit 
1. The halo of NGC 5775 therefore has the greatest spectroscopically detected extent in emission 
of any spiral galaxy halo. The vertical emission profiles will be examined in Paper II; we merely 
point out for the current purposes that, for data averaged over 10 spatial pixels, the Ha intensity 
reaches peak values, in Slits 1 and 2, respectively, of 1.8 and 6.5 xlO~^^ erg cm~^ s^^ arcscc"^, 
falling by z = 5 kpc to 7.5 and 3.8 x 10~^^ erg cm~^ s~^ arcsec"^ (NE and SW values have been 
averaged together). 

Figure 2 shows the vertical runs of the line ratios [SII] A6716/Ha, [Nil] A6583/Ha, [SII]/[NII] 
and [OIII] A5007/Ha. Because of the 85° (Irwin 1994) inclination of the galaxy, for points within 
10" of the midplane (z ^ 1200 pc in the figures), the spatial axis refiects in-plane, highly-inclined 
disk structure rather than true height above the plane, and line ratio variations are more due to 
the line of sight crossing HII regions and areas between them. This fact explains why line ratio 
minima are not always at z = kpc. 

For the first two ratios, one sees, for the most part, similar behavior to many DIG halos 
previously observed, although these measurements extend to much larger heights: a general rise 
with z, with disk values of 0.2-0.3 and halo values of [NII]/Ha reaching >1 in places. In general, 
these ratios rise more slowly with z than in NGC 891. The relatively low and constant values on 
the NE side of Slit 1 arise from the bright filament, and may simply refiect the well established 
correlation between these ratios and Ha surface brightness (Wang, Heckman, & Lehnert 1998; 
Rand 1998). In common with the East side of the halo in the spectra of NGC 891 (Rand 1997, 
1998), [NII]/Ha on the NE side of Slit 1 reaches a maximum at z = 6 and subsequently falls, 
although the maximum is reached at only z = 2 kpc in NGC 891. There is a suggestion of similar 
behavior in both lines on the S W side of Slit 1 . [S II] / [N II] is reasonably constant in the halo at a 
value of about 0.7 in both slits, but shows more variation than in NGC 891, reaching extremes of 
0.4 and 1.1. 

We examine [OIII] /Ha rather than [0III]/H/3 since the former can be determined to larger 
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z. With the exception of the NE side of Sht 2, where it remains relatively constant up to 2; = 4 
kpc, this ratio shows a clear increase with z, from values of 0.05-0.4 in the disk to as high as 
1.2 in the halo. Unlike the afore-mentioned line ratios, Slit 1 does not show a significant local 
maximum in this ratio, although the final data point on the SW side suggests that a maximum 
may have been reached. The behavior is qualitatively similar to NGC 891 (Rand 1998), except 
that those measurements extended to only 2 = 2 kpc, where a value of [0III]/H/3 = 0.8 was 
reached. Assuming optically thin gas at 10,000 K, the peak [OIIIJ/Ha of 1.2 corresponds to 
[0III]/H/3 = 3.7. 

[OI]/Ha rises from around 0.02 in the disk to 0.06 aX z = 2 kpc on the SW side of Slit 2. He I 
A5876 has been detected in Slit 2, but only for four ten-pixel averages in the disk can meaningful 
values be measured. The mean of the four He I/Ha ratios is 0.039 ± 0.003 (the uncertainty is the 
dispersion of the four values). There is no significant correlation of the ratio with Ha intensity. 

These data allow halo kinematics to be probed to unprecedented heights. Figure 3 shows 
velocity centroids of the [Nil] A6583 and [SII] AA6716, 6731 lines as a function of slit position for 
data averaged over ten pixels. The heliocentric systemic velocity, as determined from HI data, is 
Vsys = 1680 km (Irwin 1994). Once again, points at z ^ 1.2 kpc represent disk emission. The 
midplane values in both slits indicate rotation speeds roughly consistent with the rotation curve 
found by Irwin (1994). The falloff towards Vgys for z ^ 1.2 kpc is expected for a differentially 
rotating disk viewed at a not quite edge-on aspect. 

Beyond this point, though, the velocities are of the halo gas. They continue to move closer to 
Vsys with increasing z, coming to within 20 km s~^ of Vsys except on the SW side of Slit 1. The 
simplest explanation of this behavior is that the rotation speed falls with z, as suggested for NGC 
891 (Rand 1997), but here becoming roughly consistent with no rotation at the largest heights. 
The exact translation from observed to rotational velocities is complicated by the unknown 
distribution of gas along the line-of-sight. This issue will be explored in Paper III. Given the 
symmetry on either side of the two slits, it is unlikely that the kinematics are dominated by the 
tidal interaction NGC 5775 is undergoing. For Slit 1, the steeper falloff on the NE side may 
reflect a receding motion of the prominent filament. In Paper III, we will attempt to constrain 
the dependence of Vrot with z and examine whether the falloff is expected for a reasonable mass 
model, or whether other hydrodynamical effects (Benjamin 2000) may be important. 



4. Discussion 

Without carrying out detailed modeling, we can already make conclusions about the possible 
source(s) of ionization and physical conditions in the DIG that may explain the emission line 

ratios and their vertical runs. As in NGC 891 and the Reynolds layer, one can already conclude 
that pure photoionization models will not be able to explain the runs of all the line ratios. Most 
problematic are the increasing values of [O III] /Ha with z, in complete contrast to photoionization 
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models. Also of concern is the relative lack of variation in [SII]/[NII]. Hence, we are led to 
consider the two possible deviations from the pure photoionization picture discussed in §1. Both 
may be able to reproduce the trends in these two line ratios, as well as [SIIJ/Ha, [NII]/Ha, and 
[OI]/Ha. In Paper II, we wiU use data on four edge-ons (NGC 891, NGC 5775, UGC 10288, 
and NGC 4302) to ascertain whether one or the other of these models provides a more adequate 
description of the data. 

Other complications include abundances and the extent to which depletions (e.g. Howk & 
Savage 1999a) are important, as well as scattering of disk light into the line of sight by extraplanar 
dust. Absorption due to dust has been shown to be prevalent in several DIG halos (Howk k, 
Savage 1999b). Modeling of NGC 891 (Ferrara et al. 1996) suggests that scattered light is a minor 
contributor to diffuse halo Ha emission (10% at z = 600 pc) and declines with z. Nevertheless, 
scattered light may be significant in some cases. 

Apart from the excitation, the kinematics of gaseous halos and the extent to which gas can 
be expelled from the disk are relevant issues not only to the nature of the disk-halo cycle, but in 
the interpretation of QSO absorption line systems. For instance, Mgll absorption is probably 
due to both halos and disks (Charlton k, Churchill 1998); consequently possible low-z analogs of 
absorbers need to be well characterized at faint levels. 



The help of the KPNO staff is greatly appreciated. We are also grateful to L. M. Haffner, R. 
Benjamin, and R. Reynolds for reading a draft of the manuscript. 
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Table 1: Summary of the Observations 

Slit Offset along Hours of Noise"^ 

Major axis'* Integration (10~^^ erg cm~^ s~^ pix~ 
1. 32" NW 6 3l 
2 20" SE 5 3.3 



"Offsets are from the HI kinematic center determined by Irwin (1994) at R.A. 14'' 53™57.6', Dec. 3° 32' 40" (2000.0). 
''Angular pixel size is 1.38 arcsec^. An intensity of 2 x lO"'^* erg cm~^ s~^ arcsec"^ corresponds to an Emission 
Measure of 1 pc cm~® 
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Fig. 1. — Slit positions overlaid on the Ha image of NGC 5775 from Collins et al. (2000). The 
cross marks the kinematic center from Irwin (1994). Contours show the disk structure. 



Fig. 2. — Dependence of line ratios on slit position for Slit 1 (filled circles) and Slit 2 (open circles). 

Line ratios plotted are (a) [SIT] A6716/Ha, (b) [Nil] A6583/Ha, (c) [SII]/[NII] and (d) [OIII] 
A5007/Ha. Although the position axis is labeled z, note that points at z ^ 1200 pc reflect in-plane 
structure more than vertical structure. 



Fig. 3. — The intensity- weighted mean heliocentric velocities of the [Nil] A6583 and [SII] 
AA6716, 6731 lines as a function of position for (a) Slit 1 and (b) Slit 2. The dashed line indicates 
the systemic velocity. 
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